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HOMBUGLAITOEB 

= parameter as defined in text (fable 1) 

B = bottom product flow rate, moles/minute 
= parameter as defined in text, (fable 1) 

D = top product flow rate, moles/minutes 
E = Murpbree vapor plate efficiency 
f = perturbation in feed rate, mole s/minute 
E = feed rate, moles/minute 

^nu “ gain for plate n composition response to a load u: 

^nx feed composition, to feed rate 

H = , holdup molc-sj Hq, is total column holdup 

= constant as defined in text (Eqn,15) 

1-^ = reduced circulation rate 

L = liquid flow rate, moles/minute 
= tray number 
E = reflux ratio 

f = time constant; is reciprocal zero, f^,f 2 ,f^ff^ are 

first, second, third and fourth time constants 
respectively, minutes 

fg « simplified model time constet, minutes 

u = load variable, f feed rate, is feed composition 

X perturbation in composition of acetone, mole fraction 

X = mole fraction of acetone in liquid phase 
y = perturbation in vapor composition of acetone, mole 

fraotion 



Y = mole fraction of acetone in vapor pha 
Greek letters 



=5 

relative volatility 


= 

constant as defined in text (fable 1) 

Subscripts 

A 

r= 

acetone 

B 

= 

bottom product 

D 

= 

top product 

F 

= 

feed 

n 

== 

tray location, 0 is condsnscr. 



N+l is roboiler 

R 


Reboiler 

¥ 

=5 

Water 

Superscripts 

i 


first steady state 

2 

= 

second steady state 



ABSTBAOI 


!Phe open loop aynanics of a pilot scale distillation 
column was studied as part of a largo project envisaging 
direct digital control for the process. Ihe variables that 
affect the steady state operation of a distillation column 
are: feed flow rate, feed composition, reflux flow rate, 
reboiler steam pressure and feed ten^jerature. The transfer 
functions for the distillation column were obtained for step 
disturbances in feed flow rate and feed composition using a 
model proposed by ¥ahl and Harriot. This was a preliminary 
and there were numerous practical problems encountered during the ; 
column operation primarily because of lack of control. Bor 
change in feed flow rate, the proposed model predicted the 
steady state behavior fairly well but the experimental results : 
during unsteady operation did not match with the prediction 
mainly because of the assumption of constant feed rate before 
and after the disturbance was violated. Bor feed composition 
disturbance, the model did not agree with the observations at \ 

all because of the large step size vhic^ resulted in significant 
change in molar feed flow rate contrary to the assumptions. 



CHAPTEE 1 


lEJEOEUCTION 

Eynamic modelling is essential to study the unsteady 
behaviour of a process. OJhe increased concern about the 
control of industrial processes has led to a pronounced interest 
in the study of their dynamics. Eistillation is one of the 
most widely used processes in chemical industry. It requires 
a good amount of energy for its operation and often stringent 
duty is assigned to it and therefore its control within close 
limits is essential for an economical and efficient running 
of the plant • 

Models for process dynamics may be obtained theoretically 
or experimentally. A theoretical analysis is based cn the 
application of the basic physical laws governing the process. 
However, these models tend to get veiy combursome computationally 
and may require extensive experimeniation for evaluation of 
the parameters involved. In a purely experimental approach, 
the process is considered as a black box, and, the output is 
related empirically to the changes in input variables. This 
gives the dynamic behavior but does not allow an insist into 
the process, as the basic physical parameters will not appear 
in the transfer function explicitly. A more realiatic model 
is usually a combination of the two where, based on theoretical 
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knowledge one chooses a model and simplifies it in the li^t 
of the knowledge of the experimental observations and finally 
a few parameters involved are evaluated empirically. 

Distillation dynamics can be expressed, quite simply even 
for multicomponent distillation in terms of state variables 
(Toijala and G-ustaffson, 197 ' 1 ) hy writing unsteady state 
material and energy balance equations. However, such a model is 
difficult to work with, in columns having many equilibrium 
stages. The study of control systems throu^ state variable 
models for such columns would require large storage space 
and computing time. A theoretical model in all its detail is 
not very suitable, for the control of a process since the 
amount of cowpatational effort is large. The complexity of 
the process model has to be reduced, also because it is likely 
to contain unnecessary information, Toijala (19T2) and Astrom 
and Eytooff (1971) used process identifical methods to obtain 
simple models for complex systems. In process modelling for 
control, the degree of sophistication of the model is a function 
of the control theory used vdiere usually, only the basic 
structure of the process is needed. Therefore, simpler models 
are more desirable. 

In an operating distillation column the variables which 
may experience random disturbances aid which need to be 
controlled are: 



3 


1) Feed flow rate 

2) Feed compo si 1:1011 

3) Reflux flow rate 

4) Vapor flow rata 

5) Feed temperature 

However, in order to extract an appropriate model, 
planned disturbances can be given* Ihese may be 

1) Step disturbance 

2) Pulse disturbance 

3) Sinusoidal disturbance 

In this investigation, the composition response of a 
multistage pilot scale distillation column was studied for 
step changes in feed flow rate and feed composition, Phis is 
a part of a continuing project envisaged to provide digital 
control for the column. While, for control, it is necessary 
to. study closed loop dynamics, the present study is preliminary, 
and only the open loop dynamic response of the column was 
investigated. 

A survey of various distillation column models is given 
in Chapter 2, along with a description of the model chosen 
for this investigation* Chapter 3 discusses the experimental 
setup and operation of the distillation column* !Ehe results 
and discussions of the expjerimental runs are given in Chapter 4. 
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Modifications in the experimental setup are suggested here 
for a smoother operation of the column. In the last Chapter, 
concluding remarks alongwith suggestions for future -work 
are included. 



CHAPTER 2 


PISTIIMTIOH OPIUM MORELS 

Process modelling from first principles usually results 
in very complicated models. On the other hand, a theoretical 
model with simplifying assumptions, and requiring some para- 
meters to be found experimentally, usually can represent a 
fairly complicated process in a simple form. This is because 
a number of small effects, combine into a few larger overall 
effects (Astrom and Eyldioff, 1971)* The correlations should be 
made such that as rany of the limiting assumptions are absorbed 
in the parameters are possible, to increase the model accuracy. 

The simplest case would be that of a one parameter model* Por 
a constant holdup and overflow rate, Robinson and Gilliland (1950) 
have shown that the response to a feed composition change in 
terms of a single time constant. Another such simple model is 
that of Armstrong and Wilkinson (1957), ^Qie accuracy in these 
models is not hi^, but they can be used for quick calculations to 
obtain order of magnitudes, Por changes in feed composition 
and reflux flow rate, the response of the reboiler and the 
condenser has been reported by Williams et al, (1965, 1972). 

^e composition response of all the stages, to all of the 
commonly encountered disturbances, has been shown by Wahl and 
Harriot (1970), Waller (1968, 1969, 1970, 1971), and Waller 
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and Jonaason (1970, 1972), In both these models the flow 
dynamics was neglected and the composition dynamics studied. 
The two principal parameters in Waller's model are obtained 
from the number of real plates, the holdups, the flows and 
the plate gains, i,e, from steady state data only. (Ehese 
parameters, T and P, are defined for the enriching section (e) 
and stripping section (s) as follows: 


= 


(M+2) 


I- = 




where 


(1a) 


(1b) 


M = number of real plates in enriching section 
N = ITumber of real plates in stripping section 
=5 feed flow rate at steady state, Moles/sec 
^eq “ (M+N+2)llp, total hold up, moles 
is hold up per plate. 



where Q is defined as 
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with 


K 

^ ^ ^ + ■’) <1° 
1° 

'^ = ^eq ;f- + I) 5° 

JJ 

M+1 


K, 




/M+F+2 


•eq 


0 ^o 

n + 4r 

o o pO R 


P 


0.5 


0.5 


(3a) 

(31>) 

(4) 


v^ere 

= gain on plate ^ ^ =0 is conflenser and 

^ = M4-1T+1 is rehoiler 

Yg = vapor flow rate in stripping section, moles/sec 
Lg = liquid flow rate in enriching section, moles/sec 

The response curves for equal holdup in plates, condenser and 
reboiler, are given in (Waller, 1978). 

The model suggested by Wahl and Harriot (1970) considera the 
differential equations for component material balances. IHie 
four points ere disturbances can enter this process are 

1) the feed composition, sj, 

2) the feed flow rate, f 

3) the boilup rate, V 

4) the reflux rate, r 

Change in the temperature of the feed corresponds to altering 
the ref Itpc and vapor rate at the feed plate, and so is a 


1 
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special case of 3) ani 4 ) applied to a part of the colum 
and therefore is not listed, separately* 

Hie dynamic response can be expressed as 


= 



-.<■> . «.) (^1 . ^.X^l 


+ r(s)( 


•Cs)' 


(5) 


where 


x^(t) = peirturbation in composition of plate n at time t 
:^i3) - feed composition perturbation in Laplace domain 
f(s) = feed rate perturbation in Laplace domain 
vCs) = vapor rate perturbation in Laplace domain 
r(s) = reflux: rate perturbation in Laplace domain. 

^n^^^ 

^tsT » TtsT » ‘vtW “rtsT transfer functions 

for each of the disturbarces* 


A component material balance in the column gives the 
following relation between the input and output quantities 
for each of the four loads (S'igure 1 )* iDhe assuaption made is 
that the column is considered as a black box, Lhe black box 
contains the total amount of material (liquid and vapor holdups) 
which is about constant, but the average composition 


which varies. 



(r+f),(Xp+Xp) 


(B+(f+r~v)), (Xg +2:^^) 

figure 1; Distillation Golum as a Black Bsx 
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( 6 ) 

(7) 

( 8 ) 

(9) 


(I'or detailed derivation see Appendix 1) 
where 

= total liquid holdup, moles 
B = feed flow rate, molea/min 
B = bottoms flow rate, mol es/min 
. D = distillate flow rate, moles/rain 
Xj, = feed composition, mole fraction 
Xo= distillate composition, mole fraction 
Xjj = bottoms composition, mole fraction 
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The interior of the coluian is described as a set of 1+2 
well mixed tanks v3iich are maintained at steady state, and are 
upset by loads introduced into this system at some point and 
withdrawn at the others. The circulating flow consists of the 
vapor stream in one direction and the liquid stream in the 



Pigure 2: Essential Interior of a Distillation Column 

(a) low holdup, hi^ circulation rate 

(b) Hi^ holdup, low circulation rate 

The parameter which describes the ability of the system 
to maintain itself at equilibrium is the ratio of the circulation 
rate to the holdup (Eigure 2). The reoiprocal of this can be 
considered as a mixing or circulation time of the column. It 
is a capacitance, divided by the flow rate, L throu^ the 
capacitance. On the other hand, the rate at •vdiich the ^stem 
reisponds to the upaets will be given by a characteristic time 
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constant for the system, T^, Is approximately equal to 
the principal time constant of the system. Because of the 
importance of these two effects, the ratio of the two is taken 
as reduced circulation rate, 


(L/Hn,) L 

Lr = ^ 

^ /rf^ \ Xln 


(VtJ 




( 10 ) 


Linearizing and solving each of equations (6) to (9) results 


in solution of as 

5 


N+1 

r «« 


G. 


nu 


m _ n=0 

■"s “ B + D. 


(11) 


RCJ 


^OTJ 


Where 

= Gain on plate n due to load change u, u is 3^ or f or v or r 

n = 0 is condenser and n * lT+1 is the re boiler 
(Bor detailed derivation of equation (11) see Appendix 1). 
Transfer Bunctions: • 

The dynamic terms in the transfer functions for the column 
may be estimated from the two parameters T^ and 11 ^ and other 
basic physical parameters of the column. The exact transfer 
functions for the column for perturbation loads were found by 
Wahl and Harriot who solved over 3©0 cases of binary distillation 
columns. They found that the time constants for all the plates 
for all types of loads are ’ identical and the differences that 
occur are because of variations in the values of zeros. 
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Except for four time constanta 3?^, T 2 t and one zero 

T^, all other time constants were found to be ne^igible, Eor 
feed, composition load, adequate representation was obtained 
in the time domain by using the principal time •onstant and 
the third largest time constant, pro^vlding the proinr zero is 
used, IThe first and the third largest time constants were 
used because the second and fourth were effectively caiKjelled 
by the zeros. For a feed rate load, however, the first and ilie 
second time constants should be used for top of the column, since 
the third and not the second time constant was effectively 
cancelled by a zero. For the bottom of the column, only 
the principal time constant was used, as the remaining time, 
constants were effectively cancelled by the zeros. Correlation 
graphs for estimating these time constants are given in 
Appendix 2 • 

The suggested transfer functions are as follows: 


Xp(s) 




^nx^^z 

(!D,j a'+1)(f3s7TT 




(123 
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For iDOttom of coltimn 


(13) 


lor top of column. 


f (s) 


- ■ ... % 

^nf 


Ifp+TT 


f(s) 


^ 

V S+1) i 




Tr^s+i7'"n!2s+W 



Vs) 

(14) 


She gains, G-^^, can "bQ calculate.i from i|ie steady state 
compositions for two values of aiay load parameter using the 
equations given in Table 1«. The d.erivationa are given in 
Appendix 1, This procedure requires the value of . the change in 
vapor composition with change in liquid composiiafon on each 
plate, 


"i 

^ ^ S ^n . 

can be obtained from equation (16) for constant relative 
volatility and constant plate efficiency. 




E OC 


(1 + ( 1) x) 


2 + (1-B)(^) (16) 


The equations in Table 1 can be used if there is a total ■' 
condenser, equiwolal overflow, a single feed stream, and top 
and bottom product streams only* These restrictions can be 
relaxed following the procedure givaa in Appendix 1* 
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TAHOE 1: CALGULATIOH CEE PLATE COMPOSITION &AINS 


n 


Step 1 • Determination of and 
n 
0 
1 

2^ • * • 


1 

1/K 


“^e\' 


1 


, , 

°n-1~ 7^°n-2 


Common term 
0 
0 




K 


T+1 






Vrvp^ ">■ 


PeVl + : 

'“-2 L”''^ 

-K! 


s -1 

%r-2 


%+2***% 



V 


n 

0 

1 

2, . . .1^ 

%+.1 

%+ 2 ***% 


Peed composition 
0 
0 
0 

• ~E 



Peed rate 
0 
0 

0 - 


. ^-1 “ % 

y 
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Table 1 contd 


Step 2; letenaination of Condenser gain Q 


ou 


P-Da 


Por feed composition load, G = 

ox ]>fBb. 




E 


E 


Por feed rate load = --y ^ 


Step 3: Eeteimination of gain on the plates 

Por all loads: G^^^^ ~ + \ ^ou* \ %. appropriate 

loads are given above. 



lEhe distillation colunai was simulated using the model 
suggested by V/ahl and Harriot. Slhe model suggested by Waller 
could also have been used, but it has a further simplifying 
assumption that the holdup on all the plates, the reboiler and 
condenser are same. Bit for this difference both the models 
are quite similar. A comparison between the two models has 
not revealed any significanfc differences concerning the amount 
of work needed or accuracy* obtained (Waller 1978). 

In both these models, the disturbances are considered 
independently. Ihey represent models for single input - sin^e 
output system, The responses to all the disturbances can be 
obtained by superposition as shown in equation (5) » This 
will give a representation of multiple input-single output i^stem*| 
lor multiple input-multiple output system the distillation 
column can be modelled hy (Meyer et al. 1979 and Ogguniake and 
Hay, 1979). 

y(a) = G(s) n(s) + 1^(3) d(s) ( 17 ) 

where 

Y(s) = vector of output variables in Laplace domain, 
such as overhead composition, side stream 
composition etc. 

U(s) = vector of input or control variables such as 

reflux flow rate, reboiler steam pressure etc. 
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d(s) =: vector of disturbance variables in Laplace 

domain such as feed composition, feed flow rate, 
fhe elements of the matrices G(s) and &^(s) can be found 1:^ 
keeping all except one, of the elements of U(s) and ^(s) zero, 
fhe effect of that variable on the output can then be studied, 
fhese models hsve time delays incorporated in them. Lhese models 
thou^ more complex are a better representation of the disti- 
llation column. But, this requires an on-line san^ile analyzer 
in order to measure the time delays accurately. Since an on-line 
sample analyzer was not available, this model could not be used 
in the present study. 

Once the transfer functions have been obtained, the model 
con be used to study control loop strategies (Beaverstock and 
Harriot, 1973 ) . lo control top product by reflux manipulation 
the beat location of the sensor would be the top product itself. 
(¥ahl and Harriot, 1971 ) as the steady state error is minimized, 
Bor step changes in load , the optimum control sensor location 
is plate 1, lor a multiple input-miiltiple output distillation 
column, various control strategies have been suggested by 
Smith ( 1957 ) » Meyer et al. ( 1979 ), Oggunaike ( 1980 ) and 
Rosenbrock ( 1974 ). 
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EXEEBIMERa?AL WORK . 


3.1 Experimental Setup ; 

The pilot saale, multistage distillation column has been 
designed as a general purpose tower which can handle common 
systems such as methanol-water, ethanol-water, wateivacetic 
acid, benzene-toluene, acetone-water etc. The design, ia baood 
on the average properties of the compounds mentioned above* It 
can also be used to distill multicomponent mixtures of some of 
the above compounds. 

It is a sieve plate column made of stainless steel, with the 
following specifications; 


i) Number of stages 

ii) Column diameter 

iii) Stage hei^t 

iv) Thickness of tray 

v) Weir length 

vi) Number of holes/ tray 

vii) Spacing 

viii) Diameter of holes 

ix) Pitch 

A thermometer pocket has been welded to each stage for 
insertion of platinum resistance themometer . The ISiermometer 


16 

22.86 cm 
22.86 cm 
0.32 cm 
16 .76 cm 
140 

Triangular 
0,635 cm 
1.52 cm. 
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is so fixed as to locate the sensor within 3 nm above the 
sieve plate in order to sense the liquid temperature on the 
tray • Also every stage has been provided with 6 mm sample 
outlet pipe, welded near the bottom flange about 25 Eun above the 
tray. Trays 4 to 12 have been provided with feed inlet pipes 
of 2,54 cm diameter. A steel skeleton framework provides support 
for the column. The bottoms accumulator has been mounted on a 
brick foundation and the stages assembled on it. 



Condenser Specifications ; 

The condenser is a 1—2 shell and tube type. It has 
been designed for a heat duly of 37 * 5 x 10 ^ cals/hr. 


Number of tubes 
O.jD, of tubes 


59 

19 mm 



20 


Length of tube 
BWG of tube 
Tube material 
Baffles material 
Shell material 


60 cm 
16 

Stainless steel 
Mild steel 
Mild steel 


The shell was electroplated inside with Nickel and Chromium^ 


Equipment Layout : 

The distillation column is located in the TJnit Operations 
Laboratory. The location of accessories is as follows: 

First Floor: 1) Reboiler 

2) Bottom product drum 
Second Floor: 1) Feed preheater 

2) Instrument panel 
Third Floor: 1) Feed drum 

2) Condenser 
5) Reflux accumulator 

The column and the reboiler have been insulated with ^ass-wool 
and lagged with magnesia. There are no controllers in the 
setup which make the model tedious and laborious. The general 
layout is given in Figimre 5, 

Measurement of Temperature : 

Temperature is measured using platinum resistance thermo- 
meters. (instrumentation Limited, Kota). The specif ications 
are given below: 
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Temperature range 

0® to 500®C 

Thermometer element 

Platinum 

Thermometer resistance 

100 at 0®G 


140 at 100 

Maximum operating pressure 

40 kg/cm^ 

Working length 

1250 mm 

Overall length 

1410 mm 

Type 

PRT 101 

The principle used here la 

that resistance of 


materials, such as platinum is very sensitive to temperature. 

For pure materials, the change in resistance is practically 
linear over a certain range. The change in resistance is 
transduced into voltage signals. The voltage signal is then 
amplified to the required range. The resistance vs temperature 
graph and the associated voltage conversion circuit are shovni 
in Appendix 3. 

Gontrol Valves : 

Four control valves, pneumatic type (instrumentation 
Limited, Kota) have been placed to control the feed flow, reflux 
flow, bottoms flow and the steam flow* The operating air 
pressure is 3 to 15 psi. 
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Flow Measurement ; 

There are four orifice plates, one in the steam line 
(1.27 cms diameter) and one each in reflux, feed and bottom 
(0.635 cms diameter) line. The pressure drop across the 
orifice plates can be measured using a differential pressure Csp) 
cell which, in turn, can be used to calculate the flow rates. 

The appropriate reflux ratio, the maximum column capacity, 
the reboiler and the condenser duly are shown in Appendix 3. 

Semi conductor type DP cells (National semiconductors) were 
used to measure the pressure drop acaross the orifice plates. 

The chip housing was filled with silicone oil to prevent the 
working fluid from coming in contact with the cell. The 
silicone oil was separated from working fluid by using a latex 
rubber sock which allowed pressure to be transmitted but did 
not allow fluid mixing. However, the sock tended to aspture 
with prolonged usage and at times flow had to be measured directly 
by collecting the flow for a fixed time. 

Composition Measurement ; 

In binary distillation it is possible to relate the 
composition of a tray to its temperature through a T—X,Y diagram 
and hence monitoring temperature of each tray is equivalent to 
monitoring its composition. However, as a check, samples of 
liquid holdup were collected periodically from sample ports 
in 30 ml. Stoppered sampling bottles and composition checked 
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by measuring the refractive index with an Abbe Hefractometer 

This process was used only to verify the 
composition at selected positions as obtained from temperature 
measurement, 

3.2 Operating Procedure ; 

Peed consisting of a known mole percent of acetone in 
water was introduced at a predetermined flow rate into the 
column* This was done after preheating the column for about 
an hour by supplying steam to "the reboiler. Peed was introduced 
at the twelfth tray. The condenser water supply was turned on, 
the moment the feed was introduced. The column was operated 
under total reflux, for approximately an hour, to enable the 
holdup in the enriching section to build up fast, after which it 
is operated at finite reflux. The bottom product value was 
opened, vdien the feed was introduced. The feed flow rate, the 
bottoms floT^ rate and the reflux flow rate were adjusted manually 
using the pneumatic control value. This was done keeping the 
steady state values (for details see Appendix 3) . 

Samples are drawn periodically from the top product, bottom 
product and trays 1, 10 and 16, The composition was analyzed 
using an Abbe Hefractometer, The temperatures on all the 
trays were also measijred periodically. The bottoms product was 
collected in the bottoms tank. Excess of -fee top product from 
the reflux accumulator also drained into the hottcm product tank. 
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Ihe mixtiiro (top and bottom products) oolleoted In the bottoms 
drum must have the same composition as the feed.. This is 
recycled periodically using a centrifugal pump. 

After steady state operation was aohleyed, l.e. viien the 
temperatures became constant, a step change was given in feed 
rate/feed composition. The temperatures and the compositions 
were again measured till fto second steady state was achieved. 

h feed supply to the column was stopped along with the ^eam 
to the* reboiler after the second steady state was achieved. 

After closing the bottom product valve, the total column holdup 
was measured manually. 



QHAPa?BR 4 


MSIJLTS MD 3)ISCITSSIQ¥S 

The model proposed by Wahl and Harriot (1971) and 

discussed earlier was used to simulate the present distillation 

column because of its simplicily. The following assumptions 
were made ; 

1» Volumetric liquid holdup on all trays is equal* 

2. Vapor holdup is negligible compared to liquid 
holdup* 

3. Equimolar overflows of vapor and liquid hold 
(McOabe Thiele assumption). 

4. Temperature measurement at any trsy can be used 
to calculate composition using T-X,Y diagram. 

5* Reflux temperature was at bubble point. This 
assumption was not strictly correct since reflux 
was Invariably cold because of lack of temperature 
control* 

6* Plate efficiency is constant (©=0.75), 

The following experimental runs were undertaken, tvro 
each for feed rate and feed composition change. The operating 
conditions for each of the runs is listed in Table 2* 



TABLE 2: OPERATING CONDITIONS 
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TABIiE 3: GAINS BOR ADI !I1HE PLAIBS 


Plato 

Number 

Run 1 

Run 3 

Run 4 

Condenser 

0,0222 

5.5847 

0.0251 

1 

0,0046 

0.9548 

0.0058 

2 

0,0021 

0.4363 

0.0026 

3 

0,0017 

0.3631 

.0,0022 

4 

0.0017 

0.3384 

0.0022 

5 

0.0017 

0.3165 

0.0024 

6 

0>0017 

0,2470 

0.0026 

7 

0,0017 

0.1449 

0.0029 

8 

0.0017 

0.2228 

C.OOAf 

9 

0.0017 

0.0113 

0.0156 

10 

0.0024 

0.0582 

0.0119 

11 

0,0025 

0.2690 

0.0039 

12 

0.0031 

0.2511 

0.0021 

13 

0,0032 

0.0284 

0.0019 

14 

0.0023 

•~0 .0464 

0.0015 

15 

0.0016 

-0.0721 

0.0010 

16 

0.0010 

-0.0907 

0,0006 

Rs boiler 

0.0004 

-0.0979 

0.0002 
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TABLE 4; TIME COITSTAOTS EOH EOUH TYPICAL THIYS 


Plato 

Number 

b 

(min) 

*2 

(min) 

"3 

(min) 

*4 

(min) 

T 

2 

(min) 

1 

47.66 

18.05 

12.70 

10.16 

1.124 

11 

47.66 

18.05 

37.15 

29.72 

37.16 

13 

47.66 

18.05 

33.73 

26.98 

-6.72 

16 

47.66 

18.05 

38.38 

30.70 

-11.20 


30 


Because of anormous experimental difficulties, to be discussed 
later, it was not possible to get more extensive data in the 
limited time available , 


Measured values of the bottoms and distillate flow rate, 
holdups and the calculated values of the gains were used to 
calculate using equation (11) . The gain for each plate was 
calculated using equations in Table 1 and are given la 
Table 3. The details of determining the reflux ratio as well 
as molar flow rates are given, in Appendix 3. The time constants 
were found by correlations given by Wahl and Harriot, These 
correlation curves are given in Appendix 2 for convenience. The 
time constants found for four typical tisys are given in Table 4* 
Equations 12 to 15 were then used to obtain transfer functions 
for various plates as obtained below. Two plates below the 
feed and two above were chosen to demonstrate the column 
behavior. Plates 1 and 11 are above the feed and 13 and 16 are 
below it. 


The transfer functions are as follows: 

Por feed rate load, f 

.0058 (1.124 a+1) 

riate 1: = (4 7 .S8 ' CT (IV.OB 5 ^0 


Plate 


.. ,0039 (37,16 s+1) . 

TT^ irW Wit ly SPT) 


^( 3 ) 0.0019 

Plate 13: = (47,66 s+TT 





Temperature 
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T{me,hrs 

FIG 5 Temperature response for two typical trays in enriching 
section for step disturbance in feed rate (Run! ) 












T{m« ,hrs 


FIG. 8 Temperature response for two typical trays In stripping 
section for step disturbance in feed compositlonfRun 4) 




Tamper ature,®C 




FI 6. 9 Temperature response for two typical troys In enriching 
section for step disturbance in feed com posit ion (Run 4) 
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Plate 16: 



Por feed composition load, 

n QKiffS 


Plate 1: 


0,9548 (1.124 34-1) 


Plate 11: 0.2690 ( 37.16 3+l! 

”• jgtaT - UYM ' ’ a+l) ' DV.Ib ' '! 


V=) 


Plato 13: 3 ^ 


0.0284 (-,6.72 s+1) 


Blate 16: 


Pigures 4 to 9 show esperimental values of temperature as 
function of time for various plates. Also plotted on the same 
figures are the temperatures expected on the respective trays 
predicted by the model. 

It is seen from these figures that the predicted values 
agree with the experimental data only as the second stea^ state 
approaches. In between the first and second steady state, 
the experimental profiles show a dip vihich is more pronounced 
near the feed tsay. The lack of fit during the unsteady state 
operation for a stop change in feed flow rate (Pigures 4 to 7) can 
be attributed to introduction of cold feed ■v^hich lowers the feed 
flow temperature reducing vapor flow upwards and its temperature 
and increasing liquid flow downwards with lower temperature. 

This effect, however is minimal on plate 1 as expected. A 
similar observation was noted by Beaverstock and Harriot (1973) 
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even for reflux controlled column. 

Eor step disturbance in feed composition, figures 8 and 9 
show the experimental temperatures to be higher than predicted 
in the stripping section. The step change of —46.6 percent is 
really tcohi^ to allow the assumptions to hold. In particular^ 
the analysis requires that because of composition change no 
significant change in molar flow rates result. However, in the 
present case the vapor flow rates changed signif i cantly (The 
vapor flow rate in stripping section changed by over 20^ and 
in the enriching section it changed b7 over 3095) Causing 
deviations. Because of this drastic change, there was flooding 
in the lower section. This resulted in the temperature of 
the liquid (with a lower acetone content) rising up. 

The location of a sensor for a single point control should 
bo on the tray where the gain is hi^. ilso, the steady state 
error should be low. Wahl and Harriot (1971) reported that the 
optimum location of the sensor for top product control, for 
step disturbance would be in plate 1 itself. , 

Experimental Difficulties and Sources of Error ; 

1 . In the existing set up gravity feed was used. However, 
due to the vapors in the column, there was a constant locking 
up of feed. This problem was partially overcome hj venting the 
vapors from time to time from the feed line. But, 12i is caused 
a disturbance in the column. Iforeover, between the two 
succsessive venting of vapors from the feed line, there was a 




59 


continuous decline of feed rate and assumption of constant feed 
rate was tox from satisfactory. A positive di^laoemcnt pun^) 
of variable capacity should be used for the constancy of feed 
flow rate. 

2* Peed Temperature ; 

In absence of a temperature controller on the feed 
preheater, it was not possible to maintain a constant feed 
temperature. Cold feed was, therefore, introduced vdiich caused 
a thermal disturbance resijOLting in abnormal behavior as 
discussed above. 

3. ReflU3i: Rato and temperature ; 

Unless steady state is reached, the rate of reflux would 
be varying making steady state still more elusive. Also if the 
reflux temperature is not at its bubble point then it will 
cause another thermal disturbance in the column. However, 
without a temperature control on the reflux or on the condenser 
it was not possible to avoid the disturbance in the present 
study, 

4. Composition Measurement : 

!Ehe compositions in the condenser, first, tenth and 
sixteenth plates and the rebeiler were periodically checked 
using a refrac tome ter. However, except for the distillate, 
the remaining samples showed erroneous readings, (Ehis is 
because of impurities present. Refractive index measurementa. 
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being sensitive to impurities are,not reliable. It is suggestea 
that an on-line gas chromatograph should be used alongwith 
appropriate arrangements for a pure vapor sample to be analyzed 
This v/ill enable the time delays to he measured accurately. 

5* I'low Measurements! 


« 


As discussed earlier, the latex rubber sock, separating 
the silicone oil filled in the IS? cell from the wo iking fluid, 
tended to rupture making the measurements of differential 
pressure across the orifice plates inaccurate. Moreover, the 
present DP cells wore found to be very sensitive to temperature 
and mechanical vibrations. 



OOHOIITSIONS MD 1ECOMMMDA35IONS 


While the- present study provided an open loop dynamic 
model for the distillation column, there were several short 
comings in the experimental program. This being the first 
investigation on the dynamics of a pilot plant scale multi- 
stage column in this laboratoiy, it was mostly e35>loratoiy in 
nature. Successive investigations should result in more exact 
dynamic s-tucLLes after the shortcomings reported in this woifc 
have been overcome. With some control on reflux flow and its 
temperature and on feed flow and its temperature and data 
logging facility on microprocessor it would be possible to 
develop better open loop as well as closed loop dynamic M>dels 
for the column. More sophisticated multiple variahle, multiple 
time delay models, which could he easily developed after 
providing some distributed control as suggested above, will 
enable designing overall control strategy for the column. 
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Appmpix 1 

qf PnoteaaT State Bouatlons fnr iin vh^ 

Loaas ana Slrnmifi pa Model T1I.A 

(a) Peoa composition load 

^he unsteady state component- balance with tbe perturbation 
is (Pigure 1) 

P(Xp+3^) = P(X^+x^) + B(Xg+Xg) + Hj (11.1.1) 

Ihe steady state component balance is 

PXp=3X^+BXj^ (11.1.2) 

substituting Eq. 11.1.2 in Eq. 11.1.1 we get 

% ~ “ ^0 (11.1.3) 

where x^^ = dx^^dt 

(b) Peed rate load 

The unsteady state component balance for feed rate load 
is 

% = PXj. + fXj, - DX^ - Ex^ + Bijj + Bxg 

+ (11.1.4) 

Substituting 11.1.2 in 11.1.4 we get 

^ ^av ~ ~ ^o (11.1.5) 

The composition will completely define the coa^osition x, 

through the column. In this case 

Xq = f^ (x^^) (11.1.6) 

^ ^2 (II..I.7) 
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EqiiSt A1*1*6 Hnd A1#1*7 wo goi? 


where 


^0 


- ^R \v 


^R= 


(A1.1.8) 

(A1,1.9) 


Equations A1.1«3 and A1.1,5 can be rearranged to give 
dx 


+ ^av = 


■vdiere 


Load 


H, 


T = 
a 


T 


(BCjj+DC^) 

Load = Fxj, or 3^f 
Cp may be evaluated as follows: 


dXc 

^0 = dx! 


) c-^) 


av 


av 


G 


ou 


where G. 


ou 


Gain of the condenser due to load u 


u = load f , or 2 ^ 

Cj^ can also bo found similarly. 


Also 


X. 


av 




u 


R+1 


Hn 




nu 


1 


(A3. 1.10) 


(A1.1.11) 


G s= Gain on plate n due to load u. 


vfcere 
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n = plate number, 0 is condenser, F+l is reboiler 
= holdup for plate n. 


A1 .2 Peri'vation of Constants in Wahl and Harriots Model i 

iThe constants and b^ are presented in Table 1 in 
Chapter 2. This derivation has been necessiated because of 
some misprints in the original paper (Wahl and Harriot, 1971) . 
for n = 0 ^n ~ ^ ~ ^ 

for n = 1 ^n “ ^^*^^1 ~ ^ 

for n = 2 



Figure; A1.1 
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A componant balance aromd tray 1 (i-igure A1.1) for perturbation 
in composition gives 


■^e ?2 - * Vo - V = 0 

The equilibrium relation ship is 




'n 


(A1.2.1) 


(A1.2.2) 


where 


7^ = - Y 


( 1 ) 


'n 


= 

'n \ 

^ = 4^' - 4^' 

(1) (2) 

\ » \ are liquic] compositions on tray n at first 

and second steady state respectively. 

(1) (2) 

^ t are vapor compositions on tray n at first and second 

steady state respectively. 

Substituting A1.2.2 in A1.2,1 and rearranging we ^t for a 
load disturbance, u 


but 


Xo 

tv) = 


ll 

U 


^ 4] t^) - 


V2 


Therefore 


(fi) = i 


(Xq/u) 

e N 


K4 + I(- x_ X 


T 


1 


'®,) (-£) - I (-£) 

^ ^ u^ e '■u 


^^2 


e l e\ 1 


r e l e^ i e 


■3 
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■^e + I 


i.e. 


similarly f or trays 3, . .. ,% we get 


•K ^6 '^—1 ^ 

" = Vi - 7 ;kj V2 


(A1.2.3) 


For tray n = 

A component balance around feed tray (Fig. A1.2) and rearrangement 


gives 





(A1.2.4) 

(A1.2.5) 


Assuming that the load is transmitted instantly in the flow 
rGsponsc, the constants for n = found as in 

enriching section. They are as follows; 
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'll 


/a + Is, , 

I. tr — ) h_ < - ° ,. b_ 




- 7^ V2 


^ _ /a Vl + Ig, ■ H, 

^ ^-1 - TC V2 


^-1 "■ % 


~x^ 

The successive substitution of tbe solution of in the 


( 11 . 2 . 6 ) 

(I1.2i7) 

( 11 . 2 . 8 ) 


subsequent stage starting from tray 1 yields 

(V) = ^ + la (V) 

X 


(11.2.9) 


where (•■~) is the condenser gains G 


ou 


The condenser gain is found as follovrs: 
a) Feed Composition load 

The overall steady state con^onent balance for the 
perturbed column is 

B(Xg+Xjj) + 3)(X^+x^) = F(Xp+3^) (11.2.10) 

Before the column is perturbed the component balance is 


+ B Xg = FXp (11.2.11) 

Substituting Eq. 11.2.11 in 11.2,10 we get 

D(l2-) = F - B ( %-) U1.2.12) 

^ u 

Solving Bq.Al.2.9 and Eq.11, 2, 12 simultaneously for x^^^ we 
get 


Xo F - B ajj 
yg = B +■ S' B" 


(11.2.13) 
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The gain for each plate ia calculatei using Eqns* A1.2.9 anfi 
At. 2. 13. 

(h) Eor feed rate load 

The overall component balance for perturbed column is 
CE+f) Xj. = I>(Xq+x^) + (B+f) (X^+Xg) (A1.2.14) 

Substituting eq.AI. 2. 11 in Eq.A1.2.14 and solving simultaneously 
with Eq.AI. 2. 12 wo get 


The condenser gain for reflux rate load and boilup loads 
can be similarly found. , 



APPENDIX 2 
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Reduced circulation rote, Lf:^ 

FIG, A.2.t Prediction of prhrKlpal time constant 



FIG. A.2.2 Prediclion of second time constant 
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APPEEDIX 3 


of Reflux Ratio Required for Pealred 
Separation of Acetone^Water I’lix'bu.re Peed: 


Eeed 

Mole fraction acetone, Xjyp = 0.15 
Mole fraction water, ^ =0,85 
Distillate ; 

Mole fraction acetone, desired 

^A,D= 0-965 

Mole fraction water, desired, 
Bottoms 


Mole fraction acetone, desired, 

^A,B=0-°2 

Mole fraction water, desired. 



%,B“ 


Eigure: A3.1 

Let 1 be the molar flow rate of Peed, D is the top product rate, 
B is the bottoms rate 
Then mass balance equations are 
(i) Overall mass balance: 

E = D + B 

(ii) Overall component balance, acetone 


• ^•?A,E " ^ ^A,D ® ^A,B 

On solving A3. 1 and A 3.2 simultaneously we get 


(A3. 1.2) 
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D = 0.138 S' and B = 0.862F 
Por a feed rate P = 35 .38 gmole/iain, 

D =4*882 gmole/min 
B = 30,513 gmole/min, 

McCabe — Thiele procedure has been used to find the suitable 


reflux ratio by trial and error (Treybal, 1980). The operating 
line for the enriching section has a slope of B/R+1. 

Por finding the operating line of the stripping section, the 
q line is first plotted. The line is found as follov/s: 
Average molecular wei^t of feed (M ) is 
= 0.15 X 58.08 + 0.85 X 18.016 
= 24.026 


Por an acetone solution in water of 0.15 mole fraction, the 
saturation temperature 


= ®.5»C 

The molar heat of solution H = -338 BJu/lb.mol = - 187.38 cal/gmol« 
at 15°C ® 

Por a feed temperature of 25°0 

Gp^ water = '' BWlb'T = 1 Oal/gms «C 

°v’. acotc,ao= °-53 
Let the base temperature be 15®C. 

Then molar enthalpy of feed is 




^P>A ^^'"^base'O 


+ H 


a 


18.016 jj85x1x(25-15Q + 58.08 |Tl5x53x(25^15^- 187.38 
= 11.507 cal/gBiole 
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Siciilarly "tho molao? eniilialpy of salrurafod liquid 

% = 1409.578 Btu/lb mole = 781.47 ral/gmole 
latent heat of vaporization of acetone at 63.5°C is 
®s,A= '7244.89 Gal/gmole 
latent heat of vaporization of water at 63.5°C 
^s,¥ “ cal/gmole 

She molar enthalpy of saturated vapor at 63.5°c 1b 


+ Z 


A,? ^s,A ^,f ^s,¥ 


= 8437.84 cal/gmole 


% “ I^f 


8437.84 - 11.507 _ ^ 

TO7.BT"-" 7BT;Tr " 


The data has been taken from (Perry, 1973) » 

Therefore slope of q line =11 

Knowing the q line, the stripping section operating line can 
be drawn. Once the operating lines are known the number of 
stages can be found. This procedure is repeated for 2 values of 

R. 

(i) R=1 

Fumbor of stages in enriching section =13 
Number of stages in stripping section = 2 
Total number of ideal stages =15 


Taking an overall efficiency of 0.75> the real trays 
needed will be = 20 

Since the total numb ey of actual stages is 16, 

R = 1 is insufficient. , 
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(il) R = 2.5 

Number of stages in enriching section = 6 
Number of stages in stripping section =2 
fotal number of ideal stages = 8 
With a stage efficiency of £.= 0.75 this is sufficient. The 
equilibrium curve and the stage determination are shown in 
figures A3. 2 and A3.3i A pseudo equilibrium curve has been 
drawn (Treybal, 1980), and the number of trays found using this. 
The temperature versus T-J is given in Rigure A3 .4. The details 
of column design are given in Saboo and Sivaram (1980) and 
Mistry (1974). 

A3. 2 Determination of Mouid and Vapor Rlow Rates in Stripping 
and Enriching Sections and the Column Capacity ; 

The procedure is given in Treybal (1980) 

Prom Appendix A3.1 

R = 2.5 

D = 0.138 P 

B ^ 0.862 P 

Por P = 35.38 gmole/min, 

D =; 4.882 gmole/min 
B = 30,513 gmole/min 
Yg = (R +1) D 

= 15.613 gmole/min 
Dq = RD = 10.730 gmole/min. 
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- 19*151 gmole/min 

Ig = P + + Yg - Vg 

= 49.648 gmole/min 

The maximimi permissible flowrate for various sections of the 
column are as follows; 

At the top; 183.8 kg/hr 
At the feed tray: 176.4 kg/hr 
At the bottom: 221,04 kg/hr 

The detailed calculations are shown in Saboo and 
Sivaram (1980), 

Determination of Roboiler Duty : 

Date of vapors generated V_ = 19,151 gmole/min 

s 

latent heat of vaporization Hj^ = 9720 cal/gmole 
Heat duty, = 9720 x 19.151 
= 186138 eal/min 
= 738,6 Btu/min 

Latent heat of steam at 30 psig = 1164.1 Btu/lb 

*733 6 

Therefore, steam required = 

= ,634 Ib/min = .288 kg/min 

Determination of Condenser Duty ; 

Por a maximum flow rate of 150 kg/hr Saboo et al, (1980) 
have shown the calculations for condenser duigr^ 



The required heat transfer coefficient U was found 
to be U = 55.97 cal/hr. cm^ °0 

The designed TJ for the condenser (Mistry, 1974), 
is 70.30 oal/hr.om2 Therefore, the oondonser whieh is 

used for a cold reflux ia auffloient aa the flow rate uaed 
(r = 35.30 gmole/mlr) la much below the maxtnium permissible 


flow rate. 
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Mole fraction acetone in liquid /vapour phase,x,y 

FIG. A.3-6 Temperature vs x,y diagram for acetone water system 

( Othmer O F. and Benanatl R.F.,Ind. Eng. ehem. 37,(1945) p.29a 




